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I Abstract 

o 

CN Using a simple operator-norm estimate we show that the solution to the second Painleve equation within the 

J^ Ablowitz-Segur family is pole-free in a well defined region of the complex plane of the independent variable. The 

^) result is illustrated with several numerical examples. 



I — p 1 Introduction and result 

(—J Since the introduction of the Painleve equations and in recent years, due to their appearance in many different 

JJ^ models of mathematical physics there has been a growing interest in the "Painleve program" of classifying solutions, 

I ^^ I connection formulae etc. (see the book [5] as a prime example of this program). 

The six Painleve equations (and all the other equations with the Painleve property) are such that the only 

^ "movable" singularities of the general solution'' are poles, whereas all the "nastier singularities" (essential singu- 

O^ larities, branchpoints) are fixed (they depend on the particular equation under scrutiny only). For example for the 

^\^ first two equations 

C^j (PI) u"{s) = 6u{sf - s , (PII) u"{s) = su{s) + 2u{sf , (1.1) 

rn 

f^ each solution is meromorphic, with s — 00 being the essential singularity. Given a special solution to a Painleve 

^^ equation, there is a certain interest in determining the location of its poles; within this circle of ideas we may 

• • mention the work of [9] where it is shown that the tritronquee solution to Painleve I has no poles in the sector 

I arg(s)| < 47r/5 and \s\ sufficiently large. Following this, and due to its relevance in the semiclassical asymptotics 

of the Nonlinear Schrodinger equation, Dubrovin conjectured [3] that -in fact- the tritronquee is pole-free in the 

Cd whole sector (not just the distal part). Asymptotic methods [10] and numerics based on Pade approximants [11] 

provide solid evidence that the aforementioned conjecture holds true. 

The study of the distribution of poles for special solutions is in general an elusive task; in very special cases, by 
using the correspondence with an appropriate Riemann-Hilbert problem, one may use certain "vanishing lemmas" 
if the corresponding problem has sufficient symmetry (typically a Schwartz-reflection symmetry of some sort, see 
for example [1, 4]). These techniques allow at best to prove that certain solutions are pole-free on the real axis, for 
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example. When venturing into the complex plane of the independent variable, the required symmetry is lost and 
those methods cannot be applied. 

Our more modest goal here is to study the pole distribution of the so-called Ablowitz-Segur family [12] of 
solution to Painleve II; 

u"{s) ^su{s) + 2u^{s), (1.2) 

u(s) ~ KAi{s) , s — >■ +c», K G C (1-3) 

A particular and very important example of this family corresponds to the Hastings-McLeod solution k = 1; 
this solution has the property [6] that it is pole-free for real values of the variable s; this property is maintained 
for |k| < 1 while for k e M \ [— 1, 1] the solution has poles on the real axis (but is still pole-free for sufficiently large 
positive s). 

The Ablowitz-Segur family corresponds to purely imaginary solutions with behaviour as in (1.3) but k G iM; 
such solutions are oscillatory as s — )■ — cx) and decay as s — ^ -f oo and are also pole-free on the real axis. 

In both the Hastings-McLeod and Ablowitz-Segur cases (and even more general solutions) a great deal of results 
are known, including connection formulas and asymptotic form of the solutions. In particular the following is known 
(which we rephrase here in a form adapted to the context of the present paper)'' 

Theorem 1.1 (See Thm. 11.1, 11.7 in [5]) For each k £ C there exists a constant R > such that the solu- 
tions of (1-2) with behavior 1.3 are pole-free in the region sector args G [— f , f ] , |s| > R- In the special case 
K = 1 (Hastings-McLeod solution) the region is arg(s) G [~f i f] U [^, ^1 , \s\ > R. 

The question arises as to how large is the region in the complex s-plane where the solution remains free of 
poles; this question is -in principle- not easy to settle as the equation (1.2) has the Painleve property, that is, 
its solutions have only poles as singularities, but their position depend on the initial data {movable singularities). 
The result we shall prove is stated here: 



Theorem 1.2 The solutions to u"{s) ~ su{s) + 2u^{s) with the behavior 




u{s) ~ KAi{s) , s — >■ -f(X), K G C 


(1.4) 


are pole free in the whole region (the fractional power is the principal one) 




3fi(s3)> |ln|K| 


(1.5) 


In particular if \k\ = 1 the region coincides with the sector args G (— 7r/3, 7r/3). 





Remark 1.1 In general the negative real axis is excluded since S2 g iM for s G — M (no matter which determination 
of the fractional power); however, if \k\ < 1 or if k £ iR it is known that there are no poles on the negative axis. 



See later Remark 2.2. We also can replace the inequality (1.5) by 3ff(s2) > | In |k|, see Remark 2.3. 



"* The results in [5] are stated for even more general solutions of the general Painleve II ODE u" = su -{■ 2u^ -\- a (our case being 
a = 0). Much more detail is contained in the cited comprehensive book: some of those results appeared first in separate papers of the 
same authors and the specific references are contained ibidem. 



2 Proof of Thm. 1.2 

The proof is based upon the following 

Proposition 2.1 The solution of (1.2) with the asymptotics (1.3) is the second logarithmic derivative of the 
Fredholm determinant of the integral operator on H := L'^{'~i+ U 7_, | dA|) with kernel 

where X±(A) denote the indicator functions of the sets 7± := M±ic, c > (with 7+ oriented from right to left and 
7_ = ^7+ oriented from left to right), as follows 

u{sf = -^Indet f Id-KX, ) (2.2) 

ds"^ \ / 

The proof is contained in [2]; it essentially is a different representation of the Airy-operator^ on L^(M, dx) with 

integral kernel 

Ai(x)Ai'(?;) - Ai'(a;)Ai(2/) i f, r c'^yW-^o.it^) 



KAi{x,y):^ ' ' '"' _ ' ' ''' = -^ dfi dX — (2.3) 



Indeed it is shown in [[2], Thm. 3.1] that (in the present notation) 



det (ldL2(^^LJ^_) - kKs) = dot Idi2([^oo)) - k^Kai 



(2.4) 

[s,oo)/ 



where the right-hand side is the Fredholm determinant that appears in the spacing distributions of the Gaussian 
Unitary Ensemble (GUE) of random matrices [13] 



-^2(5; n) ■- det Idi2([3oo)) - 't^^Ai 



and it is related to the Ablowitz-Segur family (1.2, 1.3) by 



(2.5) 

[s,oo)J 



ds2 



uis-K)"" ■.= --^\nF2{s-K). (2.6) 



The simple idea behind the proof of identity (2.4) is that of re-expressing its restriction to the semi-interval [s, cxi) 
in Fourier space. The advantage is that, while on one hand it is unclear how to extend the representation (2.5) to 
complex values of s, on the other hand the kernel Kg of Prop. 2.1 depends analytically on s. 

Remark 2.1 The determinant of the above operator is independent of the details of the contours 7±; the usual 
choice for 7-|_ would be any (smooth) contour in the upper half plane extending to infinity in such a way that 
Ije* =11 = OdAI^""). Typically this contour extends along directions arg A = ^, ^, but we can use M -|- ic, c > 
thanks to the simple computation 

^ [i{x + icf] = --icx^ -f c^ (2.7) 



^Relative to loc. cit. we have performed an overall rotation A 1— > iA so that the phase function that was 'i?(A) := ^^> x\ becomes 

the one used in the present paper 0a;(A) := -3 — h x\ (with an overall factor of —%)■ Therefore the contours that were denoted ■^l,1r 
(for "left" and "right") have been rotated to 7+ and 7_, respectively. 



Corollary 2.1 The poles of the solutions u{s) mentioned in Prop. 2.1 coincide with the zeroes of the Fredholm 
determinant det ( Id — kK^ I as a function of s. 



It follows at once that if the operator norm of kKs can be estimated above by unity in certain regions of the s-plane, 
the Fredholm determinant is therefore nonvanishing and the corresponding solution is pole-free; and this is precisely 
what we set out to accomplish below. 

We shall follow [2] and represent the operator Kg^ as a block-antidiagonal operator relative to the natural 
splitting L'^{j+ U 7_) = L^(7+) ® L^(7_), and then use the following identity of Fredholm determinants 



det 



Id 



L^ij+Uj-} 






g 1 


J" 






det [ldi2(^_) — K^Q o J^^ 



where the two operators T, Q are defined as follows; 



(-^5) (A) 



e 2 



OsW 



(g/i)(M):=et«=(^) 



G 
d^e5«=(^)5(/i) 

dA e-5e=(^)/i(A) 
2i7r 



(2.8) 

(2.9) 

(2.10) 
(2.11) 



iR+ic ^'" M - A 

The identity (2.8) holds because (as proven in [2]) both T, Q are of trace-class in l?{j^) © L^(7_). It is also clear 
that J- and Q (between the respective spaces) have the same norms by the symmetry; we shall now make a very 
rough (but sufficient) estimate of their norms, from which the (simple) proof of Thm. 1.2 shall follow at once. 
We shall consider T for definiteness, and we start by observing that it is the composition of three operators: 



• the multiplication by e 2 ^= ^'^' on L^ ( 



ic) 



ic); 



• the Cauchy operator from L^ (M + ic) to L? (I 

• the multiplication by e~5^='-^-' on L^(M -I- ic). 

The norm of the Cauchy operator 

C ■.L^{M.-ic)^ L^{M. + ic) (2.12) 

is promptly shown to be unity. To see this fact, let ji ~ x + ic and A = i/ — ic, so that the above operator can be 
written out simply as an operator of L^(K, da;), 

f{x)Ax 



{Cf){y) := — / '' ' (2.13) 

2in J^ X - y - 2ic 

Since this is a convolution operator with the function (p{x) — 2i7r(-2ic-x) ' ^^^ representation in Fourier transform^ 
is the multiplication operator by \/2Trlf{t) where the Fourier transform of the function (p{x) is easily computed with 
the aid of the residue theorem to be 



m 



1 



-2tc 



/2n 



XR+(t)- 



(2.14) 



We abuse the notation and denote the operator by the same symbol as its integral kernel. 
'^We use the convention /(i) = — ^ f„ e~**^/(z) da;. 



Since the norm of a multiplication operator is its L°° norm and the Fourier transform is unitary, we conclude that 

ll|C||| = l. 

Thus the norm of T (and similarly Q) is bounded by the L°° norm of the multiplication operator; to achieve 
the optimal estimate we must choose appropriately c so that the line M + ic passes through the saddle point of the 
phase function 9^ = h^ -\- sX. Elementary calculations give that the saddle point is A^ := ±i^/s, arg(s) S (— 7r,7r]. 
We thus need to choose cq := ^^/s = y^cos ( ^^S£l j \Ye leave to the reader the elementary verification that 



max 



e3«^W 



e 3^(«^) max 



Q-iSAp) 



e-55R(«^) , (2.15) 



3 , 



and that the maxima are taken on at the saddle points only. It follows thus that 

Ill-Fill < ||et«^(^)|U..(R_„)||et''=(^)|U^(R+,,,) =e-i^(^*) (2.16) 

A completely parallel computation shows that 

infill <e-i'*(^*) (2.17) 

and hence 

\\\K^goT\\ < |«;|2e-t^(^^) (2.18) 

Thus, if the right side of the inequality (2.18) is bounded by 1 the operator Id — k^Q o J^ is invertible and thus the 
Fredholm determinant is nonzero. The condition \K\'^e~3^^^^^ < 1 is precisely the condition of Thm. 1.2, which is 
thus proved. 

Remark 2.2 The operator Q o T is basically the Fourier transform of the kernel X[s.oo){x)Kxi{x^y); as such it is 
shown in [13] that its norm (for s G M, of course!) is strictly less than one and hence the Fredholm determinant 
(2.8) cannot vanish if \k\ < 1. Moreover it is well known that Kai on L^(M) is a positive operator and hence for 
K € iM. we also have that the determinant of Id + |Kpi^Ai cannot vanish. 

Remark 2.3 Since the operator Q o T is trace-class, it has only discrete spectrum. Tracing the inequalities about 
the norms we see that when c > the Cauchy operator has only the null-space and then has only continuous 
spectrum, because in Fourier space it is a multiplication operator by the function e^^*"^, t > (nowhere constant). 
It then follows that the spectrum of Q o T is strictly bounded by the value in (2.18) (the maximum eigenvalue 
cannot achieve the bound since it is at the boundary of the essential spectra of the various operators involved). 
Therefore we can improve the statement of Theorem 1.2 by replacing the strict inequality by >. 

Remark 2.4 The kernel of Prop. 2.1 is an "integrable kernel" in the sense of [8] and it is thus immediately 
related to Riemann-Hilbert problem for a 2 x 2 matrix; this problem turns out to be the very .standard RHP in the 
isomonodromic approach to Painleve II, see for example [7]. Thus the estimate of the kernel Kg is nothing but 
an estimate of the norm underlying singular-integral operator associated to said RHP. We could have phrased this 
paper in that language but we opted for the current presentation in the interest of brevity because it does not require 
introducing the setup of any Riemann-Hilbert problem. 
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K = 1 (Hastings-McLeod) 



Figure 1: The poles of the Ablowitz-Segur solution with various indicated values of k and the boundary of the 
pole-free region as per Thm. 1.2. The numerics have been produced using the algorithm explained in [11] based 
on Pade approximants to the solution. It also shows that the poles on the left tend to minus infinity as k tends to 
the critical value k = 1 (of course this is only a visual cue, not any proof). 



3 Comments and conclusion 

Although very simple, the result of Thru. 1.2 seems to be not fully appreciated in the literature, to the knowledge 
of the author. It is also clear from the numeric evidence in Fig. 1 that the result is almost optimal; one may 
improve it only by a more careful estimate of the operators involved, which seems to me a rather difficult task. 

We would also like to briefly comment on the special member k = 1, the so-called Hastings-McLeod solution 
[6]; again numerical evidence strongly suggests that there are no poles also in the whole sector arg(s) e [2tt/3, Att/S] 
as remarked in the introduction (see the results of [5] summarized in Thm. 1.1), this is known for sufRciently large 
\s\ but only conjectured (and observed numerically) for \s\ finite [11]. We can offer an informal argument below 
that indicates that one is not to expect an estimate bounding the spectrum of K^ within an a-priori domain that 
excludes the point 1. 

The inspection of the numerics (for example at the top right frame in Fig. 1 (k = 0.3z)), implies that the 
eigenvalues of the operator Ks{X,^) are not confined within the unit disk as s ranges in the (left) sector args £ 
[if ' If] • '^^ ^^® ^^^^ recall that the poles that are clearly visible in the left sector of the indicated piture correspond 
to at least one eigenvalue (possibly with multiplicity) taking the value -, which is larger than one in modulus. 

On the other hand for real s G M_ the operator Kg has the same spectrum of the Airy kernel because it 
is unitarily equivalent to it: the results of [13] imply that all the eigenvalues Xj{s) are (simple and) confined in 
the interval [0, 1). Therefore, by continuity, there is a region around the negative axis where all the eigenvalues 
are confined within the unit disk; on the other hand the aforementioned numerical evidence indicates that these 
eigenvalues exit the unit disk as s ranges in the left sector. In view of the above, one could expect that there are 
values of s in the left sector at which some of the eigenvalues equal one: the fact that this is not the case, for 
|s| large is implied by Thm. 1.1, and for arbitrary |s| is suggested by the numerics. In conclusion, a proof of the 
absence of poles in the whole sector cannot follow from simple estimates directly on the norm of the operator Kg 
alone and is probably much more delicate to achieve. 
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